
  #include "stdio.h"
  #include "6cjdn.c"
  main()
  { int i;
    double ar[4][4]={ {1.0,3.0,2.0,13.0},
                             {7.0,2.0,1.0,-2.0},
                             {9.0,15.0,3.0,-2.0},
                             {-2.0,-2.0,11.0,5.0}};
    double ai[4][4]={ {3.0,-2.0,1.0,6.0},
                             {-2.0,7.0,5.0,8.0},
                             {9.0,-3.0,15.0,1.0},
                             {-2.0,-2.0,7.0,6.0}};
    double br[4][2]={ {2.0,-2.0},{7.0,3.0},
                             {3.0,2.0},{9.0,1.0}};
    double bi[4][2]={ {1.0,3.0},{2.0,7.0},
                             {-2.0,9.0},{3.0,2.0}};
    if (cjdn(ar,ai,br,bi,4,2)!=0)
      for (i=0;i<=3;i++)
        printf("x(%d)=%13.6e +j %13.6e,   %13.6e +j %13.6e\n",
               i,br[i][0],bi[i][0],br[i][1],bi[i][1]);
  }

